Whole-brain multivariate hemodynamic deconvolution for functional MRI with stability selection.

Conventionally, analysis of functional MRI (fMRI) data relies on available information about the experimental paradigm to establish hypothesized models of brain activity. However, this information can be inaccurate, incomplete or unavailable in multiple scenarios such as resting-state, naturalistic paradigms or clinical conditions. In these cases, blind estimates of neuronal-related activity can be obtained with paradigm-free analysis methods such as hemodynamic deconvolution. Yet, current formulations of the hemodynamic deconvolution problem have three important limitations: (1) their efficacy strongly depends on the appropriate selection of regularization parameters, (2) being univariate, they do not take advantage of the information present across the brain, and (3) they do not provide any measure of statistical certainty associated with each detected event. Here we propose a novel approach that addresses all these limitations. Specifically, we introduce multivariate sparse paradigm free mapping (Mv-SPFM), a novel hemodynamic deconvolution algorithm that operates at the whole brain level and adds spatial information via a mixed-norm regularization term over all voxels. Additionally, Mv-SPFM employs a stability selection procedure that removes the need to select regularization parameters and also lets us obtain an estimate of the true probability of having a neuronal-related BOLD event at each voxel and time-point based on the area under the curve (AUC) of the stability paths. Besides, we present a formulation tailored for multi-echo fMRI acquisitions (MvME-SPFM), which allows us to better isolate fluctuations of BOLD origin on the basis of their linear dependence with the echo time (TE) and to assign physiologically interpretable units (i.e., changes in the apparent transverse relaxation ΔR2∗) to the resulting deconvolved events. Remarkably, we demonstrate that Mv-SPFM achieves comparable performance even when using a single-echo formulation. We demonstrate that this algorithm outperforms existing state-of-the-art deconvolution approaches, and shows higher spatial and temporal agreement with the activation maps and BOLD signals obtained with a standard model-based linear regression approach, even at the level of individual neuronal events. Furthermore, we show that by employing stability selection, the performance of the algorithm depends less on the selection of temporal and spatial regularization parameters λ and ρ. Consequently, the proposed algorithm provides more reliable estimates of neuronal-related activity, here in terms of ΔR2∗, for the study of the dynamics of brain activity when no information about the timings of the BOLD events is available. This algorithm will be made publicly available as part of the splora Python package.

Interactions between functional networks in Parkinson's disease mild cognitive impairment.

The study of mild cognitive impairment (MCI) is critical to understand the underlying processes of cognitive decline in Parkinson's disease (PD). Functional connectivity (FC) disruptions in PD-MCI patients have been observed in several networks. However, the functional and cognitive changes associated with the disruptions observed in these networks are still unclear. Using a data-driven methodology based on independent component analysis, we examined differences in FC RSNs among PD-MCI, PD cognitively normal patients (PD-CN) and healthy controls (HC) and studied their associations with cognitive and motor variables. A significant difference was found between PD-MCI vs PD-CN and HC in a FC-trait comprising sensorimotor (SMN), dorsal attention (DAN), ventral attention (VAN) and frontoparietal (FPN) networks. This FC-trait was associated with working memory, memory and the UPDRS motor scale. SMN involvement in verbal memory recall may be related with the FC-trait correlation with memory deficits. Meanwhile, working memory impairment may be reflected in the DAN, VAN and FPN interconnectivity disruptions with the SMN. Furthermore, interactions between the SMN and the DAN, VAN and FPN network reflect the intertwined decline of motor and cognitive abilities in PD-MCI. Our findings suggest that the memory impairments observed in PD-MCI are associated with reduced FC within the SMN and between SMN and attention networks.

Expansion of a frontostriatal salience network in individuals with depression.

Hundreds of neuroimaging studies spanning two decades have revealed differences in brain structure and functional connectivity in depression, but with modest effect sizes, complicating efforts to derive mechanistic pathophysiologic insights or develop biomarkers. 1 Furthermore, although depression is a fundamentally episodic condition, few neuroimaging studies have taken a longitudinal approach, which is critical for understanding cause and effect and delineating mechanisms that drive mood state transitions over time. The emerging field of precision functional mapping using densely-sampled longitudinal neuroimaging data has revealed unexpected, functionally meaningful individual differences in brain network topology in healthy individuals, 2-5 but these approaches have never been applied to individuals with depression. Here, using precision functional mapping techniques and 11 datasets comprising n=187 repeatedly sampled individuals and >21,000 minutes of fMRI data, we show that the frontostriatal salience network is expanded two-fold in most individuals with depression. This effect was replicable in multiple samples, including large-scale, group-average data (N=1,231 subjects), and caused primarily by network border shifts affecting specific functional systems, with three distinct modes of encroachment occurring in different individuals. Salience network expansion was unexpectedly stable over time, unaffected by changes in mood state, and detectable in children before the subsequent onset of depressive symptoms in adolescence. Longitudinal analyses of individuals scanned up to 62 times over 1.5 years identified connectivity changes in specific frontostriatal circuits that tracked fluctuations in specific symptom domains and predicted future anhedonia symptoms before they emerged. Together, these findings identify a stable trait-like brain network topology that may confer risk for depression and mood-state dependent connectivity changes in frontostriatal circuits that predict the emergence and remission of depressive symptoms over time.

Denoising task-correlated head motion from motor-task fMRI data with multi-echo ICA.

Motor-task functional magnetic resonance imaging (fMRI) is crucial in the study of several clinical conditions, including stroke and Parkinson's disease. However, motor-task fMRI is complicated by task-correlated head motion, which can be magnified in clinical populations and confounds motor activation results. One method that may mitigate this issue is multi-echo independent component analysis (ME-ICA), which has been shown to separate the effects of head motion from the desired BOLD signal but has not been tested in motor-task datasets with high amounts of motion. In this study, we collected an fMRI dataset from a healthy population who performed a hand grasp task with and without task-correlated amplified head motion to simulate a motor-impaired population. We analyzed these data using three models: single-echo (SE), multi-echo optimally combined (ME-OC), and ME-ICA. We compared the models' performance in mitigating the effects of head motion on the subject level and group level. On the subject level, ME-ICA better dissociated the effects of head motion from the BOLD signal and reduced noise. Both ME models led to increased t-statistics in brain motor regions. In scans with high levels of motion, ME-ICA additionally mitigated artifacts and increased stability of beta coefficient estimates, compared to SE. On the group level, all three models produced activation clusters in expected motor areas in scans with both low and high motion, indicating that group-level averaging may also sufficiently resolve motion artifacts that vary by subject. These findings demonstrate that ME-ICA is a useful tool for subject-level analysis of motor-task data with high levels of task-correlated head motion. The improvements afforded by ME-ICA are critical to improve reliability of subject-level activation maps for clinical populations in which group-level analysis may not be feasible or appropriate, for example in a chronic stroke cohort with varying stroke location and degree of tissue damage.

Accurate Bayesian segmentation of thalamic nuclei using diffusion MRI and an improved histological atlas.

The human thalamus is a highly connected brain structure, which is key for the control of numerous functions and is involved in several neurological disorders. Recently, neuroimaging studies have increasingly focused on the volume and connectivity of the specific nuclei comprising this structure, rather than looking at the thalamus as a whole. However, accurate identification of cytoarchitectonically designed histological nuclei on standard in vivo structural MRI is hampered by the lack of image contrast that can be used to distinguish nuclei from each other and from surrounding white matter tracts. While diffusion MRI may offer such contrast, it has lower resolution and lacks some boundaries visible in structural imaging. In this work, we present a Bayesian segmentation algorithm for the thalamus. This algorithm combines prior information from a probabilistic atlas with likelihood models for both structural and diffusion MRI, allowing segmentation of 25 thalamic labels per hemisphere informed by both modalities. We present an improved probabilistic atlas, incorporating thalamic nuclei identified from histology and 45 white matter tracts surrounding the thalamus identified in ultra-high gradient strength diffusion imaging. We present a family of likelihood models for diffusion tensor imaging, ensuring compatibility with the vast majority of neuroimaging datasets that include diffusion MRI data. The use of these diffusion likelihood models greatly improves identification of nuclear groups versus segmentation based solely on structural MRI. Dice comparison of 5 manually identifiable groups of nuclei to ground truth segmentations show improvements of up to 10 percentage points. Additionally, our chosen model shows a high degree of reliability, with median test-retest Dice scores above 0.85 for four out of five nuclei groups, whilst also offering improved detection of differential thalamic involvement in Alzheimer's disease (AUROC 81.98%). The probabilistic atlas and segmentation tool will be made publicly available as part of the neuroimaging package FreeSurfer (https://freesurfer.net/fswiki/ThalamicNucleiDTI).

Comparing end-tidal CO2, respiration volume per time (RVT), and average gray matter signal for mapping cerebrovascular reactivity amplitude and delay with breath-hold task BOLD fMRI.

Cerebrovascular reactivity (CVR), defined as the cerebral blood flow response to a vasoactive stimulus, is an imaging biomarker with demonstrated utility in a range of diseases and in typical development and aging processes. A robust and widely implemented method to map CVR involves using a breath-hold task during a BOLD fMRI scan. Recording end-tidal CO2 (PETCO2) changes during the breath-hold task is recommended to be used as a reference signal for modeling CVR amplitude in standard units (%BOLD/mmHg) and CVR delay in seconds. However, obtaining reliable PETCO2 recordings requires equipment and task compliance that may not be achievable in all settings. To address this challenge, we investigated two alternative reference signals to map CVR amplitude and delay in a lagged general linear model (lagged-GLM) framework: respiration volume per time (RVT) and average gray matter BOLD response (GM-BOLD). In 8 healthy adults with multiple scan sessions, we compare spatial agreement of CVR maps from RVT and GM-BOLD to those generated with PETCO2. We define a threshold to determine whether a PETCO2 recording has "sufficient" quality for CVR mapping and perform these comparisons in 16 datasets with sufficient PETCO2 and 6 datasets with insufficient PETCO2. When PETCO2 quality is sufficient, both RVT and GM-BOLD produce CVR amplitude maps that are nearly identical to those from PETCO2 (after accounting for differences in scale), with the caveat they are not in standard units to facilitate between-group comparisons. CVR delays are comparable to PETCO2 with an RVT regressor but may be underestimated with the average GM-BOLD regressor. Importantly, when PETCO2 quality is insufficient, RVT and GM-BOLD CVR recover reasonable CVR amplitude and delay maps, provided the participant attempted the breath-hold task. Therefore, our framework offers a solution for achieving high quality CVR maps in both retrospective and prospective studies where sufficient PETCO2 recordings are not available and especially in populations where obtaining reliable measurements is a known challenge (e.g., children). Our results have the potential to improve the accessibility of CVR mapping and to increase the prevalence of this promising metric of vascular health.

Lag-Optimized Blood Oxygenation Level Dependent Cerebrovascular Reactivity Estimates Derived From Breathing Task Data Have a Stronger Relationship With Baseline Cerebral Blood Flow.

Cerebrovascular reactivity (CVR), an important indicator of cerebrovascular health, is commonly studied with the Blood Oxygenation Level Dependent functional MRI (BOLD-fMRI) response to a vasoactive stimulus. Theoretical and empirical evidence suggests that baseline cerebral blood flow (CBF) modulates BOLD signal amplitude and may influence BOLD-CVR estimates. We address how acquisition and modeling choices affect the relationship between baseline cerebral blood flow (bCBF) and BOLD-CVR: whether BOLD-CVR is modeled with the inclusion of a breathing task, and whether BOLD-CVR amplitudes are optimized for hemodynamic lag effects. We assessed between-subject correlations of average GM values and within-subject spatial correlations across cortical regions. Our results suggest that a breathing task addition to a resting-state acquisition, alongside lag-optimization within BOLD-CVR modeling, can improve BOLD-CVR correlations with bCBF, both between- and within-subjects, likely because these CVR estimates are more physiologically accurate. We report positive correlations between bCBF and BOLD-CVR, both between- and within-subjects. The physiological explanation of this positive correlation is unclear; research with larger samples and tightly controlled vasoactive stimuli is needed. Insights into what drives variability in BOLD-CVR measurements and related measurements of cerebrovascular function are particularly relevant when interpreting results in populations with altered vascular and/or metabolic baselines or impaired cerebrovascular reserve.

Enhanced top-down sensorimotor processing in somatic anxiety.

Functional neuroimaging research on anxiety has traditionally focused on brain networks associated with the psychological aspects of anxiety. Here, instead, we target the somatic aspects of anxiety. Motivated by the growing appreciation that top-down cortical processing plays a crucial role in perception and action, we used resting-state functional MRI data from the Human Connectome Project and Dynamic Causal Modeling (DCM) to characterize effective connectivity among hierarchically organized regions in the exteroceptive, interoceptive, and motor cortices. In people with high (fear-related) somatic arousal, top-down effective connectivity was enhanced in all three networks: an observation that corroborates well with the phenomenology of anxiety. The anxiety-associated changes in connectivity were sufficiently reliable to predict whether a new participant has mild or severe somatic anxiety. Interestingly, the increase in top-down connections to sensorimotor cortex were not associated with fear affect scores, thus establishing the (relative) dissociation between somatic and cognitive dimensions of anxiety. Overall, enhanced top-down effective connectivity in sensorimotor cortices emerges as a promising and quantifiable candidate marker of trait somatic anxiety.

Open access dataset of task-free hemodynamic activity in 4-month-old infants during sleep using fNIRS.

Spontaneous, task-free, hemodynamic activity of the brain provides useful information about its functional organization, as it can describe how different brain regions communicate to each other. Neuroimaging studies measuring the spontaneous activity of the brain are conducted while the participants are not engaged in a particular task or receiving any external stimulation. This approach is particularly useful in developmental populations as brain activity can be measured without the need for infant compliance and the risks of data contamination due to motion artifacts. In this project we sought to i) characterize the intrinsic functional organization of the brain in 4-month-old infants and ii) investigate whether bilingualism, as a specific environmental factor, could lead to adaptations on functional brain network development at this early age. Measures of spontaneous hemodynamic activity were acquired in 4-month-old infants (n = 104) during natural sleep using functional near-infrared spectroscopy (fNIRS). Emphasis was placed on acquiring high-quality data that could lead to reproducible results and serve as a valuable resource for researchers investigating the developing functional connectome.

Neurogenetic traits outline vulnerability to cortical disruption in Parkinson's disease.

The genetic traits that underlie vulnerability to neuronal damage across specific brain circuits in Parkinson's disease (PD) remain to be elucidated. In this study, we characterized the brain topological intersection between propagating connectivity networks in controls and PD participants and gene expression patterns across the human cortex - such as the SNCA gene. We observed that brain connectivity originated from PD-related pathology epicenters in the brainstem recapitulated the anatomical distribution of alpha-synuclein histopathology in postmortem data. We also discovered that the gene set most related to cortical propagation patterns of PD-related pathology was primarily involved in microtubule cellular components. Thus, this study sheds light on new avenues for enhancing detection of PD neuronal vulnerability via an evaluation of in vivo connectivity trajectories across the human brain and successful integration of neuroimaging-genetic strategies.

A practical modification to a resting state fMRI protocol for improved characterization of cerebrovascular function.

Cerebrovascular reactivity (CVR), defined here as the Blood Oxygenation Level Dependent (BOLD) response to a CO2 pressure change, is a useful metric of cerebrovascular function. Both the amplitude and the timing (hemodynamic lag) of the CVR response can bring insight into the nature of a cerebrovascular pathology and aid in understanding noise confounds when using functional Magnetic Resonance Imaging (fMRI) to study neural activity. This research assessed a practical modification to a typical resting-state fMRI protocol, to improve the characterization of cerebrovascular function. In 9 healthy subjects, we modelled CVR and lag in three resting-state data segments, and in data segments which added a 2-3 minute breathing task to the start of a resting-state segment. Two different breathing tasks were used to induce fluctuations in arterial CO2 pressure: a breath-hold task to induce hypercapnia (CO2 increase) and a cued deep breathing task to induce hypocapnia (CO2 decrease). Our analysis produced voxel-wise estimates of the amplitude (CVR) and timing (lag) of the BOLD-fMRI response to CO2 by systematically shifting the CO2 regressor in time to optimize the model fit. This optimization inherently increases gray matter CVR values and fit statistics. The inclusion of a simple breathing task, compared to a resting-state scan only, increases the number of voxels in the brain that have a significant relationship between CO2 and BOLD-fMRI signals, and improves our confidence in the plausibility of voxel-wise CVR and hemodynamic lag estimates. We demonstrate the clinical utility and feasibility of this protocol in an incidental finding of Moyamoya disease, and explore the possibilities and challenges of using this protocol in younger populations. This hybrid protocol has direct applications for CVR mapping in both research and clinical settings and wider applications for fMRI denoising and interpretation.

The effects of multi-echo fMRI combination and rapid T2*-mapping on offline and real-time BOLD sensitivity.

A variety of strategies are used to combine multi-echo functional magnetic resonance imaging (fMRI) data, yet recent literature lacks a systematic comparison of the available options. Here we compare six different approaches derived from multi-echo data and evaluate their influences on BOLD sensitivity for offline and in particular real-time use cases: a single-echo time series (based on Echo 2), the real-time T2*-mapped time series (T2*FIT) and four combined time series (T2*-weighted, tSNR-weighted, TE-weighted, and a new combination scheme termed T2*FIT-weighted). We compare the influences of these six multi-echo derived time series on BOLD sensitivity using a healthy participant dataset (N = 28) with four task-based fMRI runs and two resting state runs. We show that the T2*FIT-weighted combination yields the largest increase in temporal signal-to-noise ratio across task and resting state runs. We demonstrate additionally for all tasks that the T2*FIT time series consistently yields the largest offline effect size measures and real-time region-of-interest based functional contrasts and temporal contrast-to-noise ratios. These improvements show the promising utility of multi-echo fMRI for studies employing real-time paradigms, while further work is advised to mitigate the decreased tSNR of the T2*FIT time series. We recommend the use and continued exploration of T2*FIT for offline task-based and real-time region-based fMRI analysis. Supporting information includes: a data repository (https://dataverse.nl/dataverse/rt-me-fmri), an interactive web-based application to explore the data (https://rt-me-fmri.herokuapp.com/), and further materials and code for reproducibility (https://github.com/jsheunis/rt-me-fMRI).

Group-level cortical functional connectivity patterns using fNIRS: assessing the effect of bilingualism in young infants.

Significance: Early monolingual versus bilingual experience induces adaptations in the development of linguistic and cognitive processes, and it modulates functional activation patterns during the first months of life. Resting-state functional connectivity (RSFC) is a convenient approach to study the functional organization of the infant brain. RSFC can be measured in infants during natural sleep, and it allows to simultaneously investigate various functional systems. Adaptations have been observed in RSFC due to a lifelong bilingual experience. Investigating whether bilingualism-induced adaptations in RSFC begin to emerge early in development has important implications for our understanding of how the infant brain's organization can be shaped by early environmental factors. Aims: We attempt to describe RSFC using functional near-infrared spectroscopy (fNIRS) and to examine whether it adapts to early monolingual versus bilingual environments. We also present an fNIRS data preprocessing and analysis pipeline that can be used to reliably characterize RSFC in development and to reduce false positives and flawed results interpretations. Methods: We measured spontaneous hemodynamic brain activity in a large cohort ( N = 99 ) of 4-month-old monolingual and bilingual infants using fNIRS. We implemented group-level approaches based on independent component analysis to examine RSFC, while providing proper control for physiological confounds and multiple comparisons. Results: At the group level, we describe the functional organization of the 4-month-old infant brain in large-scale cortical networks. Unbiased group-level comparisons revealed no differences in RSFC between monolingual and bilingual infants at this age. Conclusions: High-quality fNIRS data provide a means to reliably describe RSFC patterns in the infant brain. The proposed group-level RSFC analyses allow to assess differences in RSFC across experimental conditions. An effect of early bilingual experience in RSFC was not observed, suggesting that adaptations might only emerge during explicit linguistic tasks, or at a later point in development.

Brainhack: Developing a culture of open, inclusive, community-driven neuroscience.

Brainhack is an innovative meeting format that promotes scientific collaboration and education in an open, inclusive environment. This NeuroView describes the myriad benefits for participants and the research community and how Brainhacks complement conventional formats to augment scientific progress.

Second-language word recognition in noise: Interdependent neuromodulatory effects of semantic context and crosslinguistic interactions driven by word form similarity.

Spoken language comprehension is a fundamental component of our cognitive skills. We are quite proficient at deciphering words from the auditory input despite the fact that the speech we hear is often masked by noise such as background babble originating from talkers other than the one we are attending to. To perceive spoken language as intended, we rely on prior linguistic knowledge and context. Prior knowledge includes all sounds and words that are familiar to a listener and depends on linguistic experience. For bilinguals, the phonetic and lexical repertoire encompasses two languages, and the degree of overlap between word forms across languages affects the degree to which they influence one another during auditory word recognition. To support spoken word recognition, listeners often rely on semantic information (i.e., the words we hear are usually related in a meaningful way). Although the number of multilinguals across the globe is increasing, little is known about how crosslinguistic effects (i.e., word overlap) interact with semantic context and affect the flexible neural systems that support accurate word recognition. The current multi-echo functional magnetic resonance imaging (fMRI) study addresses this question by examining how prime-target word pair semantic relationships interact with the target word's form similarity (cognate status) to the translation equivalent in the dominant language (L1) during accurate word recognition of a non-dominant (L2) language. We tested 26 early-proficient Spanish-Basque (L1-L2) bilinguals. When L2 targets matching L1 translation-equivalent phonological word forms were preceded by unrelated semantic contexts that drive lexical competition, a flexible language control (fronto-parietal-subcortical) network was upregulated, whereas when they were preceded by related semantic contexts that reduce lexical competition, it was downregulated. We conclude that an interplay between semantic and crosslinguistic effects regulates flexible control mechanisms of speech processing to facilitate L2 word recognition, in noise.

ICA-based denoising strategies in breath-hold induced cerebrovascular reactivity mapping with multi echo BOLD fMRI.

Performing a BOLD functional MRI (fMRI) acquisition during breath-hold (BH) tasks is a non-invasive, robust method to estimate cerebrovascular reactivity (CVR). However, movement and breathing-related artefacts caused by the BH can substantially hinder CVR estimates due to their high temporal collinearity with the effect of interest, and attention has to be paid when choosing which analysis model should be applied to the data. In this study, we evaluate the performance of multiple analysis strategies based on lagged general linear models applied on multi-echo BOLD fMRI data, acquired in ten subjects performing a BH task during ten sessions, to obtain subject-specific CVR and haemodynamic lag estimates. The evaluated approaches range from conventional regression models, i.e. including drifts and motion timecourses as nuisance regressors, applied on single-echo or optimally-combined data, to more complex models including regressors obtained from multi-echo independent component analysis with different grades of orthogonalization in order to preserve the effect of interest, i.e. the CVR. We compare these models in terms of their ability to make signal intensity changes independent from motion, as well as the reliability as measured by voxelwise intraclass correlation coefficients of both CVR and lag maps over time. Our results reveal that a conservative independent component analysis model applied on the optimally-combined multi-echo fMRI signal offers the largest reduction of motion-related effects in the signal, while yielding reliable CVR amplitude and lag estimates, although a conventional regression model applied on the optimally-combined data results in similar estimates. This work demonstrates the usefulness of multi-echo based fMRI acquisitions and independent component analysis denoising for precision mapping of CVR in single subjects based on BH paradigms, fostering its potential as a clinically-viable neuroimaging tool for individual patients. It also proves that the way in which data-driven regressors should be incorporated in the analysis model is not straight-forward due to their complex interaction with the BH-induced BOLD response.

Stability-Based Sparse Paradigm Free Mapping Algorithm for Deconvolution of Functional MRI Data.

Neuronal-related activity can be estimated from functional magnetic resonance imaging (fMRI) data with no knowledge of the timings of blood oxygenation level-dependent (BOLD) events by means of deconvolution with regularized least-squares. This work proposes two improvements on the deconvolution algorithm of sparse paradigm free mapping (SPFM): a new formulation that enables the estimation of neuronal events with long, sustained activity; and the implementation of a subsampling approach based on stability selection that avoids the choice of any regularization parameter. The proposed method is evaluated on real fMRI data and compared with both the original SPFM algorithm and conventional analysis with a general linear model (GLM) that is aware of the temporal model of the neuronal-related activity. We demonstrate that the novel stability-based SPFM algorithm yields activation maps with higher resemblance to the maps obtained with GLM analyses and offers improved detection of neuronal-related events over SPFM, particularly in scenarios with low contrast-to-noise ratio.

Voxelwise optimization of hemodynamic lags to improve regional CVR estimates in breath-hold fMRI.

Cerebrovascular Reactivity (CVR), the responsiveness of blood vessels to a vasodilatory stimulus, is an important indicator of cerebrovascular health. Assessing CVR with fMRI, we can measure the change in the Blood Oxygen Level Dependent (BOLD) response induced by a change in CO2 pressure (%BOLD/mmHg). However, there exists a temporal offset between the recorded CO2 pressure and the local BOLD response, due to both measurement and physiological delays. If this offset is not corrected for, voxel-wise CVR values will not be accurate. In this paper, we propose a framework for mapping hemodynamic lag in breath-hold fMRI data. As breath-hold tasks drive task-correlated head motion artifacts in BOLD fMRI data, our framework for lag estimation fits a model that includes polynomial terms and head motion parameters, as well as a shifted variant of the CO2 regressor (±9 s in 0.3 s increments), and the hemodynamic lag at each voxel is the shift producing the maximum total model R2 within physiological constraints. This approach is evaluated in 8 subjects with multi-echo fMRI data, resulting in robust maps of hemodynamic delay that show consistent regional variation across subjects, and improved contrast-to-noise compared to methods where motion regression is ignored or performed earlier in preprocessing.Clinical Relevance- We map hemodynamic lag using breathhold fMRI, providing insight into vascular transit times and improving the regional accuracy of cerebrovascular reactivity measurements.

Sparse coupled logistic regression to estimate co-activation and modulatory influences of brain regions.

Accurate mapping of the functional interactions between remote brain areas with resting-state functional magnetic resonance imaging requires the quantification of their underlying dynamics. In conventional methodological pipelines, a spatial scale of interest is first selected and dynamic analysis then proceeds at this hypothesised level of complexity. If large-scale functional networks or states are studied, more local regional rearrangements are then not described, potentially missing important neurobiological information. Here, we propose a novel mathematical framework that jointly estimates resting-state functional networks and spatially more localised cross-regional modulations. To do so, the changes in activity of each brain region are modelled by a logistic regression including co-activation coefficients (reflective of network assignment, as they highlight simultaneous activations across areas) and causal interplays (denoting finer regional cross-talks, when one region active at timetmodulates thettot + 1 transition likelihood of another area). A two-parameterℓ1regularisation scheme is used to make these two sets of coefficients sparse: one controls overall sparsity, while the other governs the trade-off between co-activations and causal interplays, enabling to properly fit the data despite the yet unknown balance between both types of couplings. Across a range of simulation settings, we show that the framework successfully retrieves the two types of cross-regional interactions at once. Performance across noise and sample size settings was globally on par with that of other existing methods, with the potential to reveal more precise information missed by alternative approaches. Preliminary application to experimental data revealed that in the resting brain, co-activations and causal modulations co-exist with a varying balance across regions. Our methodological pipeline offers a conceptually elegant alternative for the assessment of functional brain dynamics and can be downloaded athttps://c4science.ch/source/Sparse_logistic_regression.git.

Sleep State Modulates Resting-State Functional Connectivity in Neonates.

The spontaneous cerebral activity that gives rise to resting-state networks (RSNs) has been extensively studied in infants in recent years. However, the influence of sleep state on the presence of observable RSNs has yet to be formally investigated in the infant population, despite evidence that sleep modulates resting-state functional connectivity in adults. This effect could be extremely important, as most infant neuroimaging studies rely on the neonate to remain asleep throughout data acquisition. In this study, we combine functional near-infrared spectroscopy with electroencephalography to simultaneously monitor sleep state and investigate RSNs in a cohort of healthy term born neonates. During active sleep (AS) and quiet sleep (QS) our newborn neonates show functional connectivity patterns spatially consistent with previously reported RSN structures. Our three independent functional connectivity analyses revealed stronger interhemispheric connectivity during AS than during QS. In turn, within hemisphere short-range functional connectivity seems to be enhanced during QS. These findings underline the importance of sleep state monitoring in the investigation of RSNs.